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ABSTRACT 

We present paper six of the NIRSPEC Brown Dwarf Spectroscopic Survey, 
an analysis of multi-epoch, high-resolution (R~20,000) spectra of 25 held dwarf 
systems (3 late-type M dwarfs, 16 L dwarfs, and 6 T dwarfs) taken with the 
NIRSPEC infrared spectrograph at the W. M. Keck Observatory. With a radial 
velocity precision of ~2 km s _1 , we are sensitive to brown dwarf companions in 
orbits with periods of a few years or less given a mass ratio of 0.5 or greater. We 
do not detect any spectroscopic binary brown dwarfs in the sample. Given our 
target properties, and the frequency and cadence of observations, we use a Monte 
Carlo simulation to determine the detection probability of our sample. Even with 
a null detection result, our 1 a upper limit for very low mass binary frequency is 
18%. Our targets included 7 known, wide brown dwarf binary systems. No sig¬ 
nificant radial velocity variability was measured in our multi-epoch observations 
of these systems, even for those pairs for which our data spanned a significant 
fraction of the orbital period. Specialized techniques are required to reach the 
high precisions sensitive to motion in orbits of very low-mass systems. For eight 
objects, including six T dwarfs, we present the first published high-resolution 
spectra, many with high signal to noise, that will provide valuable comparison 
data for models of brown dwarf atmospheres. 
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1. Introduction 

Star and planet formation processes both give rise to objects in the ~1 to 20 M j up range 
(e.g., Luhman et al. 2005; Bakos et ah 2009; Marois et ah 2010; Joergens et ah 2013). 
Naively, objects of higher mass are typically assumed to form primarily via the star formation 
process and objects of lower mass are assumed to form primarily from a proto-planetary 
disk. This simplification is directly testable with a variety of approaches, both theoretical 
and observational. One straightforward observational approach is the study of multiplicity: 
do brown dwarfs come in pairs as frequently and with the same binary properties as stars? 
Although trends are suggestive among stellar visual binaries, i.e. decreasing frequency with 
later spectral types (e.g., Siegler et al. 2003), the frequency of short-period brown dwarf 
binaries is relatively uncertain (Burgasser et al. 2012). 

Binary parameter space is multi-dimensional. For a given spectral type, binaries may be 
examined in terms of companion frequency, separation, mass ratio distribution, or secondary 
mass distribution. Mazeh et al. (2003), for example, showed that for main sequence M stars 
the secondary mass distribution does not conform to a standard initial mass function (IMF) 
but instead follows a relatively flat distribution for a primary sample with M~0.7 ±0.1 
M 0 and orbital period P<3,000 days. This represents one small slice of a parameter space 
that may also be studied in diverse populations: young, intermediate age, old, metal-poor, 
clustered, etc. Specific comparisons between these samples provide a wealth of diagnostics for 
understanding the similarities and differences in formation and evolution of distinct groups 
of objects. 

In a sample of 454 F6—K3 solar-type stellar systems within 25pc of the Sun, Raghavan 
et al. (2010) found a total fraction of binaries and higher order multiples of 44%Q Among 
other results, they reconfirm the existence of the brown dwarf desert, the pronounced dearth 
of brown dwarf mass (i.e. ~10—80 M j up ) companions to stars in orbits with periods less 
than a few years (e.g., Grether & Lineweaver 2006; Metchev & Hillenbrand 2009). 

Is the brown dwarf desert the result of dynamical evolution preferentially impacting 
lower mass companions (e.g., Reipurth & Clarke 2001; Armitage & Bonnell 2002) or does it 
have more to do with poorly understood barriers to the formation of tightly bound compan- 


lr Thus, although the majority of these systems may not be multiple, the majority of the stars studied 
reside in multiple systems, as previously concluded (e.g., Abt & Levy 1976; Duquennoy & Mayor 1991). 



ions of brown dwarf mass? In a radial velocity (RV) survey with a few hundred m/s precision 
of >100 stars in the Taurus star forming region (Crockett et al. 2012), no brown dwarf com¬ 
panions to 1—3 Myr old stars have been observed (Mahmud et ah 2015, in prep), indicating 
that the existence of the desert is more likely related to formation than dynamical evolution 
in origin. Among 1—10 Myr old populations, to date only one brown dwarf-brown dwarf 
short-period (<1 year) spectroscopic binary pair has been identified (Stassun et ah 2006). 
Joergens et ah (2010) found an orbital period of 5.2 years for the young Chameleon brown 
dwarf binary Cha Ho;8 and Joergens (2008) estimates a period of >12 years for the pair 
CHXR74. However, the results of Stassun et ah and Joergens are based, respectively, on a 
survey for eclipsing systems and on a relatively small sample and thus are likely incomplete. 
Joergens estimates a binary frequency among very low mass young objects of 10^g 8 %. 

Among intermediate age brown dwarf spectroscopic binaries, Basri & Martin (1999) 
found the first brown dwarf pair, PPL 15, a 5.8 day period system, in a study of the Pleiades. 
Simon et ah (2006) studied the 100 Myr old brown dwarf binary GL 569B (Zapatero Osorio 
et ah 2004), a pair with a ~2.4 year orbit, a semi-major axis of 0.89 AU. They detected 
some evidence for a third spectroscopic component in the system, yet to be confirmed. 

RV surveys among held brown dwarfs, sensitive to binaries with periods of several years 
or less and semi-major axes of a few AU, have yielded a handful of definitive detections. In a 
sample of 59 held brown dwarfs, Blake et ah (2010) found a tight binary (a<l AU) frequency 
of 2.5/ 8 ;g%. They had previously identified and measured the orbit of the ~247 day period 
system 2MASS 0320-04 (Blake et ah 2008), independently identified on the basis of spectral 
analysis by Burgasser et ah (2008), combined their RV measurements with the astrometry of 
Dahn et ah (2008) for the ~607 day period system LSR J1610-0040, and presented RV data 
on two wide substellar pairs with periods of >10 years (2MASS J15074769-1627386 and 
2MASS J07464256+2000321). Zapatero Osorio et ah (2007) measured space motions for 
over a dozen held brown dwarfs but found no spectroscopic pairs in their sample. Burgasser 
et ah (2012) presented a solution for the spectroscopic orbit of the 148 day period pair 
SDSS J000649.16-085246.3AB, in common proper motion with the very low mass star LP 
704-48, with M9 and TO components straddling the substellar limit. Although Basri & 
Reiners (2006) indicate an overall spectroscopic binary fraction for held brown dwarfs and 
very low mass stars of 11% in their RV survey of 53 targets, only three L dwarfs in their 
sample show some level of RV variability. Of these three, 2MASS J15065441+1321060 was 
subsequently shown by Blake et ah to be non-variable and 2MASS J15074769-1627386 
and 2MASS J07464256+2000321 are long-period systems identified as binaries with imaging 
observations. Other brown dwarf pairs have been identified with imaging (e.g., Loclieu et 
ah 2007), microlensing (e.g., Bennett et ah 2008), and astrometry (e.g., Sahlmann et ah 
2013). Spectral brown dwarf binaries, systems that appear single in light of existing data 



but are spectroscopically peculiar, implicating the possible presence of a companion, may be 
numerous. Bardalez Gagliuffi et al. (2014) have identified 50 candidates, although Artigau 
et al. (2009) and Radigan et al. (2013) identified cases in which the brown dwarf binary 
candidates were instead found to have heterogeneous cloud covers. Thus it is unlikely that 
all spectral binaries have <1 year orbits, but two have been confirmed as such to date (Blake 
et al. 2008; Burgasser et al. 2012). 

No short period (<100 days) held brown dwarf spectroscopic binaries are known, but 
one intermediate age and one young system with periods of just a few days were identified by 
Basri & Martin (1999) and Stassun et al. (2006), respectively. Short period systems ought 
to be the most straightforward to identify; however, RV surveys for brown dwarf multiples 
require the world’s biggest telescopes and generous time allocations, challenging to obtain, 
and are fraught with bias. Yet without such work our understanding of substcllar multiplicity 
is skewed towards the anecdotal, and astronomers’ grasp of the basis for planetary mass 
companion formation is isolated from the context of brown dwarf and stellar mass companion 
formation. 

We report here on 11 years of dynamical observations of over two dozen held brown 
dwarfs taken at high spectral resolution in the near-infrared (1R) at the Keck 11 telescope. 
The intrinsic faintness of brown dwarfs, particularly the late L and T types, presents a 
challenge to high-resolution spectroscopic observations. However, this is the only method 
by which we may derive the RV data necessary for calculating space motions, and hence 
possible moving group or cluster membership, and the telltale RV variability of a short- 
period binary. Measurements of v sin i provide a lower limit on rotational velocity, crucial 
for understanding angular momentum evolution. Ultimately, with sufficiently precise data, 
the combination of RV versus phase together with the angularly resolved orbits for the few 
year or few tens of year period systems will yield the absolute component masses of brown 
dwarfs in binaries, invaluable for furthering our understanding of brown dwarf structure and 
evolution (Konopacky et al. 2010; Dupuy et al. 2010). Because brown dwarfs emit the 
bulk of their energy at wavelengths greater than ~1 /jrn, 1R spectroscopy provides the best 
approach for their RV measurements. 

Our goals were to identify brown dwarfs in spectroscopic binary systems and to measure 
the dynamical properties of any such pairs discovered. This project to identify short-period 
brown dwarf multiples is the latest contribution to the NIRSPEC Brown Dwarf Spectroscopic 
Survey (BDSS; McLean et al. 2001, 2003, 2007; McGovern et al. 2004; Rice et al. 2010) 
and leverages over a decade of observations to characterize brown dwarfs at high spectral 
resolution. We find that a critical factor in a productive survey hinges on the RV precision; 
sensitivity to RV variability scales rapidly with this parameter. We describe our sample, 
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observations, and data reduction in §2 and discuss our data analysis in §3. Section 4 provides 
a discussion of our results and we briefly summarize our work in §5. 


2. Observations and Data Reduction 

Targets were selected for a range of spectral types across the span of late M, L, and 
T dwarfs and on the basis of magnitude (J < 15 mag) and accessibility from the Keck 
Observatory (6 > —30°). The complete target list and observing log appears in Table 1 
which lists the object name (column 1), Right Ascension and Declination (columns 2 and 3), 
spectral type (column 4), 2 MASS J magnitude (column 5), reference for discovery paper 
(column 6), and the UT dates of observation (column 7). 

Observations were carried out with the high-resolution, cross-dispersed echclle mode of 
the facility, near-infrared, cryogenic spectrograph NIRSPEC (McLean et ah 1998, 2000) on 
the Keck II 10 m telescope on Mauna Kea, Hawaii. The NIRSPEC science detector is a 
1024 x 1024 pixel ALADDIN InSb array; a 256 x 256 pixel HgCdTe array in a slit viewing 
camera was used for source acquisition. The N3 (J-band) filter with the 12 x 07432 (3-pixel) 
slit, an echellc angle of 63.00°, and a grating angle of 34.08° produces a resolving power of 
R = A/AA ps 20,000 and nearly continuous coverage from 1.165—1.324/un (orders 58—65; 
McLean et al. 2007). Observations made on 2000 July 25 and 29 employed the 12 x 07576 
(4-pixel) slit, yielding a resolution of ~15,000. Internal white-light spectra, dark frames, 
and arc lamp spectra were obtained for flat-fielding, dark current correction, and wavelength 
calibration. Science exposures were made in 600 s nodded AB pairs at two locations along 
the slit. 

All spectroscopic reductions were made using the REDSPEC package, software produced 
at UCLA by S. Kim, L. Prato, and I. McLean specifically for the analysis of NIRSPEC 
data[^] as described in McLean et al. (2007). Wavelength solutions were determined using 
the OH night sky emission lines in each order; 4—5 OH lines across each of the orders 
used yielded wavelength solutions with typical uncertainties of better than 0.4 km s -1 . The 
two spectral orders most favorable for the analysis, 62 for the L dwarfs (1.221 pm-1.239 
/jrn; Figure 1) and 59 for the T dwarfs (1.283 fim— 1.302 pm; Figure 2), were selected 
independently on the basis of the presence of deep inherent lines in the brown dwarf targets. 
Furthermore, an additional advantage of these particular orders is the absence of terrestrial 
absorption lines, thus avoiding the necessity of division by a featureless telluric standard 
star. This provided the optimal approach for several reasons: (1) eliminating division by 


2 See: http://www2.keck.hawaii.edu/inst/nirspec/redspec.html 
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telluric standards maximized the signal-to-noise ratio and avoided the possible introduction 
of slightly offset spectra and potential small shifts in the brown dwarf absorption lines and 
hence RV measurements, (2) focusing on the narrowest and deepest lines available yielded 
the highest possible RV precision; although the K1 lines in orders 61 and 65 are deep (e.g., 
McLean et al. 2007), their breadth is unfavorable to precision RV measurements through 
cross-correlation, and (3) selecting orders 62 and 59 further guaranteed the best possible RV 
precision given the regular spacing of the OH night sky emission lines across both of these 
orders, required for a superior dispersion solution; this condition was not met for all orders 
in our J band setting. Multiple-epoch sequences for the L2 dwarf Kelu-1 and the peculiar 
T6p dwarf 2M0937 are shown in Figures 3 and 4, respectively. 


3. Analysis 
3.1. Radial Velocities 

Rice et al. (2010) found typical systematic RV uncertainties of 1—2 km s -1 for a sample 
observed with a similar methodology, similar signal to noise, and with some overlap in target 
data with this paper (Table 2). We thus adopt a conservative 2 km s -1 internal uncertainty 
for our RVs. We tested this estimate by cross-correlation of the RV invariant target 2M0036 
(Blake et al. 2010), an L4 dwarf, for which we obtained 7 epochs over more than 5 years. 
The maximum RV shift between epochs was 1.91 km s -1 ; the standard deviation in the RV 
shift for all epochs was 0.59 km s -1 . Thus 2 km s -1 provides a reasonable if conservative 
internal uncertainty on individual RV measurements. 

At least two, and as many as seven, spectra were taken for each of our targets. We 
tested for radial velocity variability by cross-correlating the highest signal-to-noise spectrum 
against the spectra from all other epochs for a given target; no significant variability was 
detected (§4). Table 2 lists the number of spectra taken for each object (column 2) and the 
total number of days spanned by the observations (column 3). RVs (column 4) were either 
taken from Blake et al. (2010) or determined by cross-correlation of the highest signal-to- 
noise spectrum for a particular object with spectra of objects with known RV (from Blake 
et al.) and similar spectral type, sometimes of type both earlier and later than our target. 
The RVs resulting from cross-correlation of a target with more than one other object were 
averaged and the standard deviation added in quadrature with the internal uncertainties in 
the radial velocity measurements. This result in most cases was dominated by the 2 km s^ 1 
internal uncertainty; however, for a few objects, primarily the fainter and thus lower signal 
to noise late T dwarfs, this procedure resulted in an uncertainty of 3 km s^ 1 (Table 2). 
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We use the average RV values from Blake et al. (2010) when available because of their 
unprecedented precision, obtained by fitting models to near-IR K-band CO bandhead at ~2.3 
/jrn target spectra. The models are composed of synthetic template spectra pins observed 
telluric spectra; the CO bandhead region of the telluric spectrum is rich in deep CH 4 lines 
that provide a wavelength dispersion and zero-point reference with a precision as good as 
a few tens of m/s. Small iterations of the RV shift between the synthetic photospheric 
spectra and the telluric spectra allow for high accuracy in the target RV measurements. We 
compared our results with values from other RV studies in the literature, e.g., Basri et al. 
(2000), Mohanty & Basri (2003), and Zapatero-Osorio et al. (2007). In every case our RVs 
were comparable to other values within 1 a. We provide our results where indicated in Table 
2 . 


3.2. Rotational Velocities 

Column 5 of Table 2 lists the v sin i values for our targets. Most of these were taken 
from the literature (Basri et al. 2000; Mohanty & Basri 2003; Rice et al. 2010; Blake et al. 
2010; references given in column 6 ). To estimate v sin i values for the remaining targets, we 
used visual comparison with objects of neighboring spectral types after superimposing the 
spectra. For some objects we convolved comparison spectra of known v sin i with a boxcar 
kernel in order to produce resulting spectra of larger v sin 1 for comparison. This method was 
approximate and yielded uncertainties of 5—10 km s -1 , based on visual comparisons with 
objects of known v sin i, for the T dwarfs in our sample. Nevertheless, these are the first 
estimates available for some of the targets and thus provide a useful guide. 


4. Discussion 

4.1. Field Brown Dwarf Spectroscopic Multiplicity 

Konopacky et al. (2010) obtained angularly resolved imaging and spectroscopy for 
each component in 24 very low mass stellar and brown dwarf subarcsecond visual binaries 
contributing to eventual measurements of orbital solutions and component masses. Our goal 
was to use high spectral resolution observations to identify any RV variability over time 
that might indicate a spectroscopic binary brown dwarf. This requires binaries with orbital 
periods sufficiently short to measure the component motion at a significant level, i.e. at least 
several a above the RV uncertainty. 

To explore our sensitivity to the brown dwarf binary parameter space, given our ~2 km s ^ 1 



RV precision, we ran a Monte Carlo simulation of 10 5 possible binary orbits for each of the 25 
objects in our sample, following Burgasser et ah (2014). Orbital parameters were uniformly 
distributed in log semi-major axis (10 -3 to 10 2 AU), mass ratio (0.8—1.0), sine inclination 
(0—1), eccentricity (0—0.6; Dupuy & Liu 2011), and all other orbital parameters (argument 
of periapsis, longitude of ascending node, and mean anomaly, 0—27r). We converted our 
target spectral types to effective temperature using the empirical relation of Stephens et al. 
(2009), and then from effective temperature to mass using the evolutionary models of Bur¬ 
rows et ah (2001), assuming ages of 0.5 Gyr, 1.0 Gyr, and 5 Gyr. Each of the 10 5 simulated 
orbits was sampled at the dates given and the primary orbital RV was calculated. A binary 
detection, for a given semi-major axis bin (0.2 dex), was defined as a system for which a 
maximum RV difference between all dates was >3cr, i.e. >6 km s^ 1 given our 2 km s^ 1 
precision. The results are summarized in Figure 5. The most important factor impacting 
the probability of detecting a potential binary was the frequency of observation for a given 
target (Table 2). 

A binary with a separation of <0.1 AU should in principle be straightforward to detect 
with an RV precision of 2 km s _1 . However, given our estimated target masses and sampling 
frequency, and assuming an age of 1.0 Gyr, we could have detected such an orbit only 50% 
of the time for only 12 of the sources in our sample (middle left-hand panel of Figure 5). 
The detection probability for an 0.1 AU binary fails to reach 90% for any of our sources. 
Using the probabilities of detection for separations greater than a given threshold a, P(> a), 
as a measure of the effective sample size, N e ff(a) = S,P(> a), we find our null result 
translates into a la upper limit of 18% for spectroscopic binaries down to a = 0.1 AU, based 
on binomial statistics. Only for systems with separations below 0.01 AU (~1 day orbits) 
could the spectroscopic binary frequency of our sample be characterized as relatively rare, 
i.e. <10%. 

These limits apply when we consider the detectability of individual systems. However, 
a signature of unresolved multiplicity could also emerge in higher velocity dispersions for 
the sample as a whole. Identifying higher dispersions across the sample requires robust 
determination of the individual measurement uncertainties, but we can perform a rough 
assessment as follows. Using the same simulation parameters, we calculated the distribution 
of velocity dispersions one would obtain if a given fraction of sources (randomly selected) 
were binaries with semi-major axes in logarithmically spaced bins. For a sample devoid of 
binaries, the mean dispersion is somewhat less than the adopted measurement uncertainty, 
about 1.75 km s _1 . Sources with radial orbital motion drive the mean velocity dispersions 
of the sample higher. Figure 6 displays the thresholds at which the mean simulated sample 
velocity dispersions are 1.5, 3 and 5 times higher than the dispersions assuming a 2 km s _1 
measurement uncertainty. The most conservative threshold is reached at a semi-major axis 
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of 0.03-0.04 AU, and is detectable at even small binary fractions (i.e., 1-2 sources in the 
sample being binary). This analysis is roughly consistent with the individual detection limits 
above, and again implies that we can rule out a significant fraction of binaries (> 10 %) only 
for separations <0.01 AU. 


4.2. Notes on Known Visual Binaries 

Of the 25 targets in our sample, 7 are known visual binaries. For these systems we 
estimated the upper limit for the observable RV shift, A(RV) maa; , for the brighter binary 
component between two epochs, assuming the most favorable possible observing conditions: 
( 1 ) the epochs correspond to the two phases at which the primary component is moving 
toward and away from us with maximum RV, (2) the projected separation corresponds to 
the semi-major axis of the system, and (3) the orbit is circular and edge-on. The observed 
and estimated binary properties are given in Table 3. A discussion of each visual binary and 
the results of our observations follows. 


4.2.1. 2MASS J22344161+4041387 - M6 

Using laser guide star adaptive optics imaging at the Keck II telescope, Allcrs et al. 
(2009) identified 2M2234 as a 1 Myr year old, visual binary with a projected physical sepa¬ 
ration of 51 AU. Given the observed binary properties, the A(RV) max is ~1.9 km s _1 (Table 
3). With an orbital period of 824 ^ 3(0 years the inclination is effectively indeterminable. 
This estimate for the period is based on a circular orbit. A more realistic value, IOOO^oq 0 
years, is calculated in Allers et al. (2009). In either case, it is not possible to observe the 
system at phases separated by half the orbit. Furthermore, given the v sin i of 17 km s ^ 1 
for 2M2234 (Table 2), it is also impossible to spectroscopically resolve the RVs of the two 
components in a single epoch spectrum, even though the component near-IR magnitudes are 
almost equal, because the maximum relative component RV separation is significantly less 
than the rotational broadening. 

Cross-correlation of our 3 epochs of spectra with each other demonstrated no RV shift 
between the 2007 and 2009 data. Between the 2006 and 2007 data there was an apparent shift 
of —7.7 km s -1 ; however, the signal to noise of the 2006 spectrum (~20) is considerably lower 
than that of the other epochs (~80), and the spectra are veiled (Allers et al. 2009), thus we do 
not have confidence in the 2006 result. Using the young M 6 brown dwarf [GY92] 5 for cross- 
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correlation with our 2007 spectrum, we obtain an RV of — 10±2 km s _1 (Table 2)[^J similar to 
the results of Allcrs et al. on the basis of Keck HIRES data from 200(|^J — 10.6±0.5 km s _1 . 
Allers et al. raise the possibility that 2M2234 could be a higher order multiple system, 
which would account for the overluminous nature of the A component. Our multi-epoch 
observations failed to detect any short-period, i.e. P <a few years, hierarchical spectroscopic 
binary in this system, although our sensitivity to intermediate separation binaries, and binary 
orientations unfavorable for detection, limit any significant statistical conclusions (§4.1). 
Given the greater K s — L' excess in the 2M2234A, it is feasible that the excess luminosity 
is related to the circumstellar disk structure, orientation, and/or possible accretion activity. 
Such a mismatch in disk properties around the components of very low mass binaries is not 
unprecedented; for example, the TWA 30AB wide, co-moving pair has an apparently edge-on 
disk around the embedded, earlier-type component, extinguishing this late type M star by 5 
magnitudes with respect to the cooler component (Looper et al. 2010). 


4.2.2. 2MASS J07464256+2000321 - LI 

The 2M0746 binary is a nearby (d ~ 12 pc), tight (~3 AU) system. We use the 
Konopacky et al. (2010) astrometric measurements (Table 3) to determine a A(RV) max of 
2.0 km s -1 . Konopacky et al. find an average primary/secondary flux ratio of 1.5 ±0.1, 
challenging the assumption that angularly unresolved spectra are fully dominated by the 
primary (Blake et al. 2010). 

We observed 2M0746 at two epochs separated by almost exactly 4 years, about 1/3 of 
the orbital period. Cross-correlation of our two order 62 J-band spectra yielded a 1.3 km s^ 1 
shift with a high correlation coefficient, 0.92. Comparing the epochs of our observations with 
the RV curve plotted for this system in Figure 14 of Blake et al. (2010), this is almost exactly 
the expected result; however, we are not sufficiently confident in our RV uncertainties to give 
it much weight. 


3 Cross-correlating the same 2M2234 spectrum with another young M6, CFHT Tau 7, Rice et al. (2010) 
found —13.1 km s -1 . 

4 This result is the weighted mean of two RVs; Shkolnik et al. (2012) use the same Keck measurements 
to determine an unweighted mean of —10.9±0.7 km s _1 . 
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4.2.3. Kelu-1 - L2 

A rapid rotator with v sin i of ~70 km s _1 and Li absorption, Kelu-1 was identified as 
a brown dwarf by Ruiz et al. (1997). Martin et al. (1999) hypothesized that Kelu-1’s over¬ 
luminosity and Li abundance might be explained by a young age or an additional component 
in the system (e.g., Golimowski et al. 2004). Liu & Leggett (2005) using Keck AO imaging 
found that Kelu-1 was a 0 / . / 291 binary. Gclino et al. (2006) estimated spectral types for 
the components of L2 and L3.5 and a total mass of 0.115 ±0.014 M 0 . In an unpublished 
preprint, Stumpf et al. (2008) describe additional observations of the system with VLT AO 
imaging through 2008; the separation steadily increased to 07366 in 2008. The position angle 
has not changed by more than 4° or 5°. Adopting the largest separation observed by Gclino 
et al. as the semi-major axis, 07298 ±07003, we estimate a period of 39 ±5 years based on 
a circular orbit (although Stumpf et al. favor a high eccentricity of 0.82 ±0.10). If viewed 
edge-on, this implies a A(RV) maa , of 4.3 ±0.4 km s -1 , marginally detectable with our ~2 
km s -1 precision. 

Measurements of the Kelu-1 system RV in the literature are inconsistent: Basri et al. 
(2000) found 17 ±1 km s _1 in June of 1997 and Blake et al. (2010) determined RVs of 
6.35 ±0.39 and 6.41 ±0.75 km s -1 in March and April of 2003. On the basis of angularly 
resolved spectra of the two known components, Stumpf et al. (2008) suggest that Kelu-1 A 
is itself a spectroscopic binary. We used our highest signal to noise (S/N) ratio spectrum 
of Kelu-1 to cross-correlate against five other epochs (Figure 3), all of S/N ratio >50 per 
resolution element (the January, 2006, spectrum, with a S/N of ~10, was not included in this 
analysis). Our RV measurements, from 2002 through 2011, show RV shifts of <3 km s -1 . 
We did not detect any clear evidence in our spectra for additional motion resulting from the 
A-component moving in a relatively short-period spectroscopic orbit; however, this could 
conceivably be the result of binary properties and/or viewing geometry (§4.1). 


4.2.4. 2MASS J15074769-1627386 - L5.5 

Over a 6-year baseline, Blake et al. (2010) detect a marginally significant (<2 0 ) trend in 
the RV of 2M1507, a nearby (d=7.3 pc) L dwarf. They obtain a false alarm probability of 2.2 
% and suggest the possibility that 2M1507 is a >5000 day binary with an angular separation 
of 074. However, deep, high-resolution imaging sensitive to a contrast ratio of 5 magnitudes 
(Bouy et al. 2003; Reid et al. 2006) has not revealed any companions. No significant RV 
variations are evident in the 5 high-resolution spectra we obtained between 2000 and 2008; 
cross-correlation of the highest S/N ratio spectrum (LIT 2000 April 25) against the other 4 
epochs resulted in RV shifts of <1.7 km s -1 with an uncertainty of ~2 km s -1 . This result, 
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however, does not rule out multiplicity; Blake et al. observed ~0.5 km s -1 of motion over 
6.5 years, thus we would not expect much more than that over our 8 year baseline. Given 
the lack of definitive evidence for multiplicity, this system is not included in Table 3. 


12.5. DENIS-P J0205.4-1159 - L5.5 

Koerner et al. (1999) initially identified this system as binary. Bouy et al. (2005) 
describe evidence for a third object in a bound orbit with the secondary component. The 
estimated properties of the wide binary orbit are uncertain but the period is at least 47 years 
and the A(RV) max is at most 4.4 km s^ 1 (Table 3). For the presumed close binary, Bouy et al. 
estimate an orbital period of 8 years and a semi-major axis of 1.9 AU, implying a A(RV) maa , 
of ~7 km s -1 . Our three spectra of DENIS 0205, taken in 2001 and in 2006, are of low S/N 
ratio. Cross-correlation between the epochs yields —2.7 and —2.1 km s -1 with a correlation 
coefficient of only ~0.4, reflecting the poor quality of the data. Sufficiently frequent and 
deep imaging and RV monitoring of this system may provide the requisite phase coverage, 
preferably with better precision than 2 km s -1 , to determine a full orbital solution for the 
inner binary over the course of one orbital period. 


4.2.6. DENIS-P J1228.2-1547 - L6 

Using the Hubble Space Telescope , Martin et al. (1999) identified DENIS 1228 as the 
first angularly resolved brown dwarf - brown dwarf pair with a separation of 0 / . / 275±0'. / 002 
(Bouy et al. 2003). After several years of monitoring the components’ positions, Brandner 
et al. (2004) estimated the orbital properties of the system, listed in Table 3. The A(RV) maa , 
for this binary, 4.3 km s -1 , in combination with the period of ~44 years from Brandner et al., 
is not favorable for the detection of an RV shift over the 4 year time scale of our NIRSPEC 
observations. Cross-correlating our 2007 May spectrum with those taken in 2011 February 
and June yields a 0 km s -1 RV shift. Continued monitoring of the visual orbit with high 
angular resolution imaging and high precision RV spectroscopic techniques will help to refine 
the parameters over the next decades, necessary to determine individual component masses 
in the long term. 
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4-2.7. SDSSp J042348.57-041403.5 - TO 

SDSS 0423 is one of the visual brown dwarf binary systems which spans the L and 
T classes. Burgasser et al. (2005) measured a separation of 0 / . / 16 and estimated a total 
mass of 0.08—0.14 M 0 . Assuming that the separation is equal to the semi-major axis of the 
system, 2.50±0.07 AU, the period falls in the range of 10.5 to 13.9 years (Table 3) and the 
A(RV) maa; is 5.3—7.1 km s -1 . We observed the system in 2001, 2005, and 2006, covering 
close to half of the estimated orbital period. However, the cross-correlation between the 2001 
and 2005 spectra yielded a shift of only —0.4 km s -1 and between 2001 and 2006 of 1.77 
km s -1 , indistinguishable within the uncertainty of our RV measurements, especially because 
the 2006 spectrum was particularly noisy. Thus we find no evidence for significant orbital 
motion, implying a longer period or an unfavorable viewing geometry for the detection of an 
RV shift, or both. 


4.3. 2MASS J05591914-1404488 

The T4.5 dwarf 2M0559 presents an enigmatic case of an over-luminous, extremely low- 
mass object. Observers and theorists alike have speculated (Burgasser 2001; Dahn et al. 
2002; Burrows et al. 2006; Dupuy & Liu 2012) that this system is an equal mass binary. 
Alternatively, there may be fundamental processes at play in the mid-T dwarf atmospheres 
that are not yet well-understood. Specifically, this source is the lynchpin in the J-band 
brightening/cloud disruption scenario (Burgasser et al. 2002a). Zapatero Osorio et al. 
(2007) estimate limits on possible planetary mass companions in this system, but such a 
secondary component would not explain the unusually high brightness. 

We obtained four observations of 2M0559 over a time baseline of 6.5 years. For an 
age of 1 Gyr, our Monte Carlo simulation (§4.1) indicates a 50% detection probability for 
a threshold semi-major axis of 0.13 ALI and a 90% detection probability for a threshold 
semi-major axis of 0.003 AU. The threshold semi-major axis is the separation below which 
a spectroscopic companion would be detected with a particular probability. Burgasser et al. 
(2003) rule out the presence of a relatively bright companion object closer than 0709. At the 
~10 pc distance to 2M0559 (Dahn et al. 2002), 0709 corresponds to ~0.9 ALL Thus, ample 
parameter space for a bright binary companion to this object remains unexplored and our 
confidence in a null result for a companion object is only high (>90%) for extremely short 
periods of days or less. Monitoring this system with extremely precise RV measurements (see 
next section) with regular cadence over a considerable time baseline will help to fill in the 
potential binary parameter space gap and might also provide insight into the atmospheric 
properties. 



4.4. The Importance of High Precision RV Measurements 


For spectroscopic binary systems, Figure 7 illustrates the relationships between the 
primary object’s mass, the primary orbital velocity, and the orbital period on the basis 
of Kepler’s third law. We show results for three distinct values of the mass ratio (q); a 
circular, edge-on orbit is assumed for simplicity. For a system with a primary of mass 0.08, 
the substellar limit, and a mass ratio of 1.0, the primary object’s RV is ~3.5 km s” 1 for a 
period of 12 years, approximately the shortest period system among the visual binaries in 
our sample (Table 3). With a precision of 2 km s -1 , motion of the primary (or the secondary, 
given a mass ratio of unity) in such a system is only detectable for very specific phases and 
viewing angles. The probability of detection with 2 km s^ 1 precision increases for shorter- 
period binaries; however, again, this is only true under certain specialized conditions (§4.1). 
None of the multi-epoch spectra in our sample of 25 brown dwarf systems reveals more than 
~3 km s^ 1 of RV variability. Even for the seven known brown dwarf binaries observed, some 
with a cadence that regularly sampled a significant fraction of the estimated orbital period, 
we were unable to unambiguously detect any RV variability. 

Specialized techniques for the highly precise measurement of small RV shifts, such as 
those applied to high-resolution K-band spectra by Blake et al. (2007, 2010), Prato et al. 
(2008), Konopacky et ah (2010), Bailey et ah (2012), Burgasser et ah (2012), and others, are 
required to reliably detect motion in brown dwarf binaries, even for those with orbital periods 
as short as days. In their 6-year study of late-type M and L dwarfs with NIRSPEC on the 
Keck II telescope, Blake et ah (2010) obtained a precision of 200 m s -1 on slowly rotating 
L dwarfs, providing sensitivity to orbital motion of brown dwarf binaries with periods of 
decades and mass ratios as low as ~10% (Figure 7), the upper limit for the detection of 
giant planetary companions. 

In the study described here, even for our sample of 25 systems with zero detections 
of spectroscopic binaries, it is still not possible to use the results to definitively character¬ 
ize short-period low mass binaries as rare. The sampling and geometry of such systems 
are simply not well-suited to identification with our 2 km s -1 precision and random ob¬ 
serving cadence. Thus as far as it is possible to say with the extant data, very low mass 
spectroscopic binaries are not necessarily intrinsically rare, but even with one of the largest 
samples available, statistics show (§4.1) that 2 km s^ 1 uncertainties provide relatively weak 
constraints. 
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5. Summary 

We obtained multiple-epoch spectra of a sample of 25 very low-mass field dwarfs, three 
M dwarfs, sixteen L dwarfs, and six T dwarfs, between 2000 April and 2011 June to search 
for RV variability and spectral evidence for multiple components. With a precision of ~2 
km s -1 , we were sensitive to RV variability at a statistically significant level only in systems 
with periods of about a day or less, assuming a favorable distribution of orbital properties 
and viewing geometries relative to our line of sight. In none of the systems studied, including 
the seven known, wide binaries observed, did we detect any RV variability >3 km s -1 . For 
over a dozen objects in our sample we present the first published high-resolution spectra 
and provide RVs and rotational velocities for the entire sample, either based on this work 
or taken from the more precise measurements in Blake et al. (2010). We show multi-epoch 
spectral sequences for two objects of particular interest, Kelu-1 and 2M0937, an L2 and 
a peculiar T6, respectively. No significant variations are seen in these or the other target 
spectra, some of which boast an exquisite S/N ratio in excess of 100. 

RV measurements of brown dwarfs are important both for the ultimate measurement of 
brown dwarf masses (Konopacky et al. 2010) and for the spectroscopic detection of very low- 
mass, even planetary, companions to presumed single brown dwarfs (Blake et al. 2010). The 
close binary fraction of very low mass systems is highly uncertain (e.g., Bardalez Gagliufk et 
al. 2014). We conclude with the observation that to satisfy these scientific goals requires high 
S/N ratio, strategic sampling cadence, and relatively high precision measurements: with the 
200 m s -1 precision of Blake et ah, it is possible to detect several-Jupiter mass companions 
even in orbits of decades (bottom panel, Figure 7). Long-term monitoring programs of binary 
brown dwarfs, and in particular candidate spectroscopic binary brown dwarfs (Bardalez 
Gagliufk et ah), with high spectral resolution, component-resolved spectroscopy (Konopacky 
et ah), with high spectral resolution unresolved spectroscopy (Burgasser et al. 2012), and 
with high-angular resolution imaging (e.g., Radigan et al. 2013), over time scales of days to 
years are required. Results of these efforts will yield component mass measurements with 
sufRcient precision to stringently test models of brown dwarf structure and evolution, and, 
in the case of younger systems, formation (e.g., Schaefer et al. 2014). It is crucial that 
RV monitoring programs take advantage of high-precision techniques for a future high-yield 
science return. 
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Fig. 1.— The NIRSPEC echcllc order 62 M and L dwarf spectral sequence for BDSS RV 
targets. Spectra are corrected for barycentric velocity, normalized, boxcar smoothed to a 
3-pixel resolution element, and offset by a constant. 
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Fig. 2.— The NIRSPEC echellc order 59 T dwarf spectral sequence for BDSS RV targets. 
Spectra are corrected for barycentric velocity, normalized, boxcar smoothed to a 3-pixcl 
resolution element, and offset by a constant. 
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Fig. 3.— Multiple-epoch spectra for NIRSPEC echclle order 62 of the visual binary Kelu-1. 
Spectra are corrected and shifted as in Figures 1 and 2. We did not measure any significant 
RV shift over the nearly nine year baseline of the observations (§4.1.3). 
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Fig. 4.— Multiple-epoch spectra for NIRSPEC echclle order 59 of the peculiar T6 dwarf 
2MASS 0937+2931. Spectra are corrected and shifted as in Figures 1 and 2. No significant 
RV shift or line variability was observed in these spectra. 
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Fig. 5.— Results from our Monte Carlo simulation showing the number of our targets 
with 50% (left panels) and 90% (right panels) spectroscopic binary detection probability 
as a function of the threshold semi-major axis (§4.3), given our target properties and RV 
precision. Because mass is a function of age, we present results for 3 age bins; spectral 
types were converted to mass estimates using models and empirical relations (see text). The 
outliers Kelu-1 and 2M0036, more favorable to detection because of the combination of their 
relatively high mass and our high sampling frequency, are not included in the right hand 
panels but ranged between 0.011—0.015 AU for a 90% detection probability. Kelu-1 is not 
included in the left panels, either, but its detection probabilities for the three ages ranged 
from 0.34—0.47 AU for a 50% detection probability. In addition to Kelu-1 and 2M0036, only 
15/25 objects had 90% detection probabilities, thus fewer objects are plotted in the right 
panels. Bin sizes are 0.025 AU on the left and 0.0008 AU on the right. 
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Fig. 6.— Binary fraction detection thresholds as a function of the typical semi-major axis 
based on the velocity dispersion of our sample. Colored contours correspond to velocity 
dispersions 1.5, 3, and 5 times higher (from right to left) than the observed dispersion 
cr, computed by sampling randomly oriented orbits with the given semi-major axis at the 
observational epochs of our sample. The most conservative threshold (5cx) corresponds to 
semi-major axes of 0.03-0.04 AU. 
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Fig. 7.— Each panel illustrates the relationship in a spectroscopic binary between primary 
object mass (Ma), primary orbital velocity (iq), and orbital period, curved lines, for 3 sample 
mass ratios, q=l, 0.5, and 0.15, top to bottom. Plots for periods of 50 years (solid line), 20 
years (short dash line), 12 years (dash-dot line), 5 years (long dash line), and 1 year (dotted 
line) are shown. The horizontal dotted line demarcates the nominal sub-stellar cut off for 
the more massive star in the system. The orbits are assumed to be circular and edge-on, the 
most favorable possible orientation for detection. 
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Table 2. Radial and Rotational Velocities 


Object 

(1) 

N 

Obs. 

(2) 

AT 

(days) 

(3) 

<RV> 

(km s -1 ) 

(4) 

v sini 
(km s^ 1 ) 

(5) 

References 
< RV >, a sini 
(6) 

GJ 406 

3 

1488 

194=2 

<84=2 

Rice et al. (2010) 

2M2234+40 a 

3 

2166 

—104=2 

174=2 

This work, Rice et al. (2010) 

2M0140+27 

3 

566 

94=2 

114=2 

Rice et al. (2010) 

2M0345+25 

3 

1864 

64=3 

304=5 

This work 

2M0746+20 a 

2 

1470 

52.374=0.06 

32.724=0.56 

Blake et al. (2010) 

2M0208+25 

3 

3056 

204=2 

124=2 

This work 

2M1658+70 

2 

499 

—25.604=0.12 

12.264=0.76 

Blake et al. (2010) 

2M2057—02 

3 

499 

—24.684=0.43 

60.564=0.37 

Blake et al. (2010) 

2M0015+35 

2 

3054 

—37.354=0.16 

10.234=2.55 

Blake et al. (2010) 

Kelu-l a 

7 

3214 

6.374=0.35 

68.884=2.60 

Blake et al. (2010) 

2M2104—10 

2 

893 

—21.094=0.12 

23.444=0.23 

Blake et al. (2010) 

G196-3B 

4 

2159 

—24=2 

104=2 

This work, Mohanty & Basri (2003) 

2M1506+13 

4 

796 

—0.684=0.11 

11.394=0.94 

Blake et al. (2010) 

GD165B 

5 

3214 

—294=2 

184=2 

This work, Mohanty & Basri (2003) 

2M2224—01 

2 

2500 

—37.554=0.09 

25.494=0.41 

Blake et al. (2010) 

2M0036+18 

7 

1965 

19.024=0.15 

35.124=0.57 

Blake et al. (2010) 

2M1507—16 

5 

2887 

—39.854=0.05 

21.274=1.86 

Blake et al. (2010) 

DENIS 0205 a 

3 

1554 

74=2 

224=5 

This work, Basri et al. (2000) 

DENIS 1228 a 

3 

1347 

24=2 

224=3 

This work, Basri et al. (2000) 

SDSS 0423—04 a 

3 

1554 

284=2 

604=10 

This work 

SDSS 1254-01 

2 

1478 

44=3 

154=5 

This work 

2M2254+31 

2 

709 

144=3 

154=5 

This work 

2M0559—14 

4 

2353 

—94=3 

104=5 

This work 

2M2356—15 

2 

144 

194=3 

154=5 

This work 

2M0937+29 

5 

1487 

—54=3 

604=10 

This work 


a Binary system. 



Table 3. Observed and Estimated Binary Properties 
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